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Abstract 

We study the BTW-height model of self-organized criticality on a 
square lattice with some long range connections giving to the lattice the 
character of small world network. We find that as function of the fraction 
p of long ranged bonds the power law of the avalanche size and lifetime 
distribution changes following a crossover scaling law with crossover ex- 
ponents 2/3 and 1 for size and lifetime respectively. 

Small world networks have recently been observed to exist in many physi- 
cal, biological and social cases ||l|, H, ^, ^ |^. A simple minded example is the 
structure of the neo-cortex and the network of acquaintances in certain soci- 
eties. Another ubiquitous phenomenon (SOC) ^, ^ in nature and society 
is self-organized criticality, i.e. the appearance of avalanches of all sizes with- 
out characteristic scale over a certain range. Easily one can imagine cases in 
which both of these phenomena occur simultaneously, that means that we find 
an avalanche type spreading on a sparsely long-range connected network. As 
an example we propose the spreading of neural information inside the cervical 
cortex which is due to the threshold behaviour given by the firing rule which au- 
tomatically induces avalanches of synapses. Another example one can imagine 
are societies where each individual has a certain threshold of endurance. Let us 
take the example of Peter Grassberger for the classical SOC-model of office 
clerks moving sheets of papers from desk to desk, as illustrated in Peter Bak's 
book in Fig. 13 and admit that the clerks do not sit on a square lattice but 
have a more realistic connectivity in their work relation (typically small world 
network behaviour) . Many other examples of this kind can be thought of. 

In this short paper we want to present a model for self-organized criticality on 
a graph having the properties of small world networks. In fact, we consider the 
classical height model of Bak, Tang and Wiesenfeld (BTW) ||^ on a square lattice 
which has been "rewired" for a certain fraction of bonds p which are chosen of 
arbitrary range. For that purpose we take a square lattice of linear size L and 
all bonds present between nearest neighbours sites. Then we choose randomly 
two sites of the system and place a bond between them (which can therefore be 
of any length smaller than ^/2L). In order to keep the coordination of the sites 
on average to be four, one of the smaller bonds going to a neighbouring site of 
one of the end points of one long bond is removed. This procedure repeats until 
a fraction p of all bonds has been replaced by long range connections. This type 
of graph is not the same as the one used by most authors working on small world 
networks, because the underlying short range lattice is not a linear chain but 
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a square lattice. Nevertheless the properties should qualitatively be the same: 
large world behaviour for large p and a small world behaviour for small p. The 
situation p = corresponds to the simple square lattice and p = 1 to random 
graph with average coordination 4 and long range connections (Viana-Bray) on 
which one would expect mean-field behaviour. 

On each site of the lattice we place an integer value hi less than a threshold 
he = At each time step randomly one site i is chosen and its value hi is 
increased by unity, i.e. a unit mass is added to it. 

When the value of the height hi having q neighbours reaches the threshold 
/ic, it topples, i.e. it distributes mass equally to its neighbours: 

hi hi - q (1) 

hj — !■ hj + I 

where j goes over all q neighbors of site i. We see that eq. (1) preserves the 
mass (the sum of all heights). 

We have studied the statistics of avalanches monitoring as well their size s 
(number of sites that toppled at least once) as their lifetime t, that means the 
number of time step an avalanche lives. The quantities n(s) and n(t) denote 
the avalanche size and life time distributions. For the case p = 0, that means 
the classical BTW-model on the square lattice it is known that asymptotically 



nis) ^ ,s-^ (2) 
n{t) - r'' 

with a — 1.0 and r = 0.5 in two dimensions. 

It is our aim here to investigate what happens for different values of p. To 
this purpose we have analysed square lattices of size L — 50, 100 and 200 and 
a range of p values from p = to p = 0.5. The lattice has spiral boundary 
conditions in one direction and two open boundaries at top and bottom, where 
mass in excess can flow out of the system. In each configuration we have injected 
randomly 2500 particles of unit mass, letting each time the avalanche proceed 
until no site had a value oi h > he- The data were averaged over 200, 50 and 
10 configurations for L — 50, 100 and 200 respectively. 

In Fig. [| we see the distribution of avalanche sizes for different system sizes 
and two different values of p, 0.1 and 0.3 in a double logarithmic plot. We see 
that the data follow a straight line for nearly two decades indicating that we 
still find SOC-behaviour. The slope of the straight line gives us the exponent 
a. In Fig. |^ we see the corresponding figures for the distribution of life times. 
Again the data show power law behaviour and the slope gives the exponent r. 
We observe that the exponents a and r depend on p but they do not appear to 
depend on L. 

In Fig. ^ we see the dependence of a and r on p as obtained from our 
simulations. We see that for p — + we obtain the classical result on the square 
lattice of BTW ct = 1, r = 1/2 §. For p close to unity the values of the 
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exponents convert to the meanfield values cr = 3/2 and r = 1 |T^, This 
is of course not surprising and the question is if the continuous change of the 
exponent a and r is an intrinsic continuous hne of critical points or, if we have 
here a crossover phenomenon as it appears for instance in magnetic models that 
interpolate between 2 and 3 dimensions. For that purpose we tried for both 
the size and the lifetime avalanche data in Fig. ^, a data collapse of all the 
distributions for different values of p following the classical crossover scaling 

ripis) - s-''T{sp*) (3) 

np{t) = t-^g{tp^) 

where and G are scaling functions and and ip are universal crossover ex- 
ponents. We see from Fig. ^ that a collapse of the data works reasonably well 
yielding crossover exponents of = 0.65 ±0.1 for the avalanche size distribution 
and -0 = 1.0 ± 0.1 for the lifetime distribution. 

By studying the BTW-model on a small world square network we observed 
a crossover to meanfield behaviour following a crossover scaling law of eq. (3). 
It would be interesting to see if the same occurs for the simpler Manna-model 
and we have heard calculations are already under way JT^ . 

This work has been partially supported by the European TMR Network- 
Fractals under contract No.FMRXCT98G183 and by MURST-PRIN-2000. 
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Figure 1: Log- log plot of the distribution of avalanche size n(s) as function of s 
for different system size L and for p = 0.01(a) and p = 0.3(5). The continuous 
lines indicate the slopes 0.98 and 1.45 for (a) and (b) respectively. The statistics 
is of 2500 particles injected in 200,50 and 10 configurations for L — 50, 100 and 
200 respectively. 




Figure 2: Log-log plot of the distribution of avalanche lifetime n{t) as function 
of t for different system size L and for p = 0.01(a) and p — 0.3(6). The con- 
tinuous lines indicate the slopes 0.54 and 0.99 for (a) and (6) respectively. The 
statistics is the same as in Fig. |l|. 



Figure 3: Critical exponents a and r for the distribution of avalanehe size 
and lifetime as function of p. The exponent are calculated for the system size 
L = 50. 
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Figure 4: Data collapse of the distributions of avalanche size n{s) and lifetime 
n{t) for different values oip and L ~ 50. (a) The plot of the quantity n{s)s'^ for 
a = 1.0 as function oisp'^ gives a value for the crossover exponent (j) = 0.65±0.1. 
(b) The analogous plot of the quantity n{t)t'^ for r = 0.5 as function oftp^ gives 
a value for the crossover exponent ijj = 1.0 ± 0.1. 
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